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Abstract. Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. We developed a 
model-independent method for determining the WIMP mass by using data (i.e., measured recoil energies) of direct detection 
experiments. Our method is independent of the as yet unknown WIMP density near the Earth, of the form of the WIMP 
velocity distribution, as well as of the WIMP-nucleus cross section. It requires however positive signals from at least two 
detectors with different target nuclei. At the first phase of this work we found a systematic deviation of the reconstructed 
WIMP mass from the real one for heavy WIMPs. Now we improved this method so that this deviation can be strongly reduced 
for even very high WIMP mass. The statistical error of the reconstructed mass has also been reduced. In a background-free 
environment, a WIMP mass of ~ 50 GeV could in principle be determined with an error of ~ 35% with only 2 x 50 events. 

Keywords: Dark Matter, direct detection, WIMP mass 
PACS: 95.35.+d, 29.85.Fj 



INTRODUCTION 



srf is defined as 



By now there is strong evidence that more than 80% of 
all matter in the Universe is dark (i.e., interacts at most 
very weakly with electromagnetic radiation and ordinary 
matter). The dominant component of this cosmological 
Dark Matter must be due to some yet to be discovered, 
non-baryonic particles. Weakly Interacting Massive Par- 
ticles (WIMPs) x with masses roughly between 10 GeV 
and a few TeV are one of the leading candidates for Dark 
Matter (for reviews, see Refs. [1]). 

Currently, the most promising method to detect many 
different WIMP candidates is the direct detection of the 
recoil energy deposited in a low-background laboratory 
detector by elastic scattering of ambient WIMPs on the 
target nuclei [2]. The differential rate for elastic WIMP- 
nucleus scattering is given by [1]: 
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where po is the WIMP density near the Earth and Ob is 
the total cross section ignoring the form factor suppres- 
sion. The reduced mass m r .N is defined by 
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where m x is the WIMP mass and mjn that of the target 
nucleus. Finally, v m j n = CCy/Q is the minimal incoming 
velocity of incident WIMPs that can deposit the energy 
Q in the detector and 
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Here R is the direct detection event rate, i.e., the number 
of events per unit time and unit mass of detector material, 
Q is the energy deposited in the detector, F(Q) is the 
elastic nuclear form factor, f\ (v) is the one-dimensional 
velocity distribution function of the WIMPs impinging 
on the detector, v is the absolute value of the WIMP 
velocity in the laboratory frame. The constant coefficient 



In our earlier work [3], we developed methods for re- 
constructing the one-dimensional velocity distribution, 
/l(v), and for estimating its moments from the recoil 
spectrum as well as from measured recoil energies di- 
rectly in direct detection experiments: 
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Here r thre = (dR/dQ)Q = Q tja<t is an estimated value of the 
scattering spectrum at the threshold energy, gthre, and /„ 
can be estimated by 
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where the sum runs over all events in the data set. Note 
that Eq.(5) can be used without prior knowledge of the 
local WIMP density, p , of the velocity distribution func- 
tion of incident WIMPs, f\ (v), as well as of the WIMP- 
nucleus cross section, Cq. 



DETERMINING THE WIMP MASS 

By requiring that the values of a given moment of f\ (v) 
estimated by Eq.(5) from two detectors with different tar- 
get nuclei, X and Y, agree, we found a simple expression 
for determining the WIMP mass [4]: 
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Here n^O, m (xj) an d F(x,y)(Q) are me masses and 
the form factors of the nuclei X and Y, respectively, 
and f"thre,(x,F) refer to the counting rates at the threshold 
energies of two detectors. Note that the form factors in 
the estimates of I„ x and I nY using Eq.(6) are different. 

Additionally, for the spin-independent (SI) scattering 
of a supersymmetric neutralino, which is the perhaps 
best motivated WIMP candidate [1], and for all WIMPs 
which interact primarily through Higgs exchange, the SI 
scattering cross section is approximately the same for 
both protons p and neutrons n. Writing the "pointlike" 
cross section oq of Eq.(2) as 

°b = (|)m r 2 N A 2 |/ p | 2 , (9) 

where f p is the effective XXPP four-point coupling, A is 
the atomic number of the target nucleus, one finds 
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The factor S appearing in the denominator is the expo- 
sure of the experiment, which is dimensionless in natural 



units; it relates the actual counting rate to the normalized 
rate (1). Since the unknown factor po|/p| 2 on the left- 
hand side here is identical for different targets, Eq.(10) 
leads to another expression for determining m x . 
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Here we have assumed m^ X Y ) A^ XY y an d introduced 
the quantity 
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In order to yield the best-fit value of m x as well as its 
statistical error, we introduced a % 2 function: 

X ? E {fix -Ay) % l (fj,x-fjj)- (13) 

It combines the estimators for different n in Eq.(7) with 
each other, and with the estimator in Eq.(l 1). In Eq.(13) 
we defined 
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we analogously defined also n max +2 functions fij. Here 
"max determines the highest moment of f\ (v) that is in- 
cluded in the fit. The are normalized such that they are 
dimensionless and very roughly of order unity; this alle- 
viates numerical problems associated with the inversion 
of their covariance matrix. The first « max + 1 functions 
are basically our estimators of the moments in Eq.(5). 
The last function is essentially the right-hand side in 
Eq.(10). It is important to note that m x in Eqs.(14a) and 
(14b) is now a fit parameter, not the true (input) value of 
the WIMP mass. Recall also that our estimator (6) for /„ 
appearing in Eqs.(14a) and (14b) is independent of m x . 
Hence the first n max + 1 fit functions depend on m x only 
through the overall factor a 1 . 

Moreover, ^ here is the total covariance matrix. Since 
the X and Y quantities are statistically completely inde- 
pendent, ^ can be written as a sum of two terms: 
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The entries of this matrix involving only the moments of 
the WIMP velocity distribution can be read off Eq.(82) 
of Ref. [3], with an obvious modification due to the nor- 
malization factor in Eq.(14a). Since the last defined in 
Eq.(14b) can be computed from the same basic quanti- 
ties, i.e., the counting rates at gthre and the integrals Iq, 
the entries of the covariance matrix involving this last fit 
function can also be computed straightforwardly. 



MATCHING THE CUT-OFF ENERGIES 



The basic requirement of our method for determining 
m x given in Eq.(7) is that, from two experiments with 
different target nuclei, the values of a given moment 
of the WIMP velocity distribution estimated by Eq.(5) 
should agree. This means that the upper cuts on f\ (v) 
in two data sets should be (approximately) equal. 1 Since 
Vcut = aVfimax. it requires that 

6max,Y = (^~^J Gmax.X- (16) 

Note that a defined in Eq.(4) is a function of the true 
WIMP mass. Thus this relation for matching optimal 
cut-off energies can be used only if m x is already known. 
One possibility to overcome this problem is to fix the 
cut-off energy of the experiment with the heavier target, 
minimize the^ 2 ( '?2^,rec) function defined in Eq.(13), and 
estimate the cut-off energy for the lighter nucleus by 
Eq.(16) algorithmically. 

In Figs. 1 we show our numerical results for the re- 
constructed WIMP mass based on Monte Carlo sim- 
ulations. 76 Ge and 28 Si have been chosen as two tar- 
gets. The dotted (green) lines show results for Qmax.Ge = 
Gmax.Si = 100 keV with a systematic deviation discussed 
in Ref. [4], whereas the solid (black) lines have been ob- 
tained by using Eq.(16) with 2m.ax.Ge = 100 keV and the 
true (input) WIMP mass m x ; n . The dashed (red) lines are 
for the case that 2 ma x.Ge = 100 keV, whereas 2 ma x,Si has 
been chosen such that X 2 ( m x,rec) I s minimal. As shown 
here, with only 50 events (upper frame) from one ex- 
periment, our algorithmic process seems already to work 
pretty well for WIMP masses up to ~ 500 GeV. 

SUMMARY 

In this paper we described the basic ideas of our method 
for determining the WIMP mass and gave the main for- 
mulae. The algorithmic process for correcting the sys- 
tematic deviation of the reconstructed WIMP mass has 
also been discussed. More details and discussions about 
determining m x can be found in Ref. [6]. 
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Here we assume the threshold energies to be negligibly small. 



50 + 50 events. Si and Ge, standard halo. C- < 100 keV 
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FIGURE 1. Results for the reconstructed WIMP mass as 
well as its error interval based on the combined fit, using 
Eq.(13) with « max = 2. We assume that the scattering cross sec- 
tion is dominated by spin-independent interactions, The the- 
oretical predicted recoil spectrum for the shifted Maxwellian 
velocity distribution function [1], [3] with Woods-Saxon elastic 
form factor [5], [1] (vo = 220 km/s, v e = 231 km/s) have been 
used. We simulated with 50 (upper) and 500 (lower) events on 
average from each experiment before cuts. 
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